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Abstract. The dynamical percolation transition of two dimensional Axial Next 
Nearest Neighbour Ising (ANNNI) model to pulsed magnetic field has been studied 
by finite size scaling analysis (by Monte Carlo simulation) for various values of 
frustration parameters, pulse width and temperature (below the corresponding 
static transition temperature). It has been found that the size of the largest 
geometrical cluster shows a transition for a critical field amplitude. Although the 
transition points shift, the critical exponents remain invariant for a wide range of 
frustration parameters. It is also same as that obtained for the 2d Ising model. 
This suggests that although the static phase diagrams of these two models differ 
significantly in various aspects, the dynamical percolation transition of both these 
models belong to the same universality class. 



PACS numbers: 64.60.ah, 05.70.Fh, 61.43.Bn 



2 



1. Introduction 

The response of magnetic systems to time dependent external magnetic fields has been 
of current interest in statistical physics ([l]-[T2]). Magnetisation reversal dynamics is 
an important subject for many technological applications of magnetic materials. In 
magnetic recording industry fast magnetisation reversal in storage media is essential 
for high data transfer rate [13] . These spin systems driven by time-dependent external 
magnetic fields, have basically got a competition between two time scales: the time 
period of the driving field and the relaxation time of the driven system, giving rise to 
interesting non-equilibrium phenomena. Magnetisation reversal by pulsed magnetic 
field has been studied extensively ([T]-[I2]). The magnetisation reversal were studied 
below the corresponding static transition temperature (T®), where the system remains 
ordered. A pulse of finite strength (h p ) and finite duration (At), applied along 
the opposite direction of the prevalent order tends to reverse the magnetisation. 
Depending on the temperature (T) and pulse duration (At), there will be a critical 
field strength (hp) at which the system will be forced to change from one ordered 
state (magnetisation too) to another ordered state (magnetisation —too). Below this 
critical strength, the magnetisation will eventually remain unchanged as the field pulse 
is taken off. This transition is called "magnetisation-reversal transition" . 

A percolation transition is also associated with this magnetisation reversal 
induced by the external magnetic field. By percolation transition we mean the 
transition of the largest "geometrical cluster" (consisting of nearest neighbor parallel 
spins) size from a large value to zero. In the initially ordered state, large clusters 
of particular spin orientation ("up" or "down") exists. But due to the pulsed field 
directed along the opposite direction, the spins gradually flip (depending on the 
magnitude and duration) and for a critical value of the field or pulse duration the 
pre-existing large clusters vanish, giving rise to oppositely oriented clusters. Thus 
if we measure the largest "geometrical cluster" size for a particular orientation at 
a certain temperature (below the critical temperature) and finite duration of field 
pulse, it shows a transition at a critical field magnitude. Although thermal percolation 
transition of the "geometrical" clusters, in 2d Ising model has been studied extensively 
[HI [T51 [TBI IT?] , very little is known about the field induced dynamical percolation 
transition [18 . It was found numerically (based on Monte Carlo simulation) that for 
pure two dimensional Ising model at finite temperature the field induced magnetisation 
transition and the percolation transition occur for the same critical magnetic field. 
Moreover the critical exponents were found (using finite size scaling analysis) different 
from that of the thermal percolation transition, indicating the dynamical percolation 
transition to belong to a separate universality class [T8,. This transition was also 
studied for the 2d Ising model with diagonal next nearest neighbor interaction for 
a typical value of the frustration parameter [18 and it was found that the critical 
exponents remain unchanged. 

Here we intend to study the field induced percolation transition for the 2d Axial 
Next Nearest Neighbour Ising (ANNNI) model which is one of the simplest classical 
Ising model with a tunable frustration and has been studied for a long time (see 
[H2 [20j EU [22] for review). The 2d ANNNI model is a square lattice Ising model with 
nearest-neighbor ferromagnetic interaction along both the axial directions and second- 
neighbor anti-ferromagnetic interaction along one axial direction. The Hamiltonian 
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Figure 1. Schematic phase diagram of the two-dimensional ANNNI model 
according recent studies 23 24 

is 

W = — J /J 8 x,y[ a x+l,y + Sx,y+1 ~ Ks x+2,y) (1) 

where the sites (x,y) run over a square lattice, the spins s x , y have binary states (±1), 
J is the nearest-neighbor interaction strength and n is a parameter of the model. 
For positive values of k the second-neighbor interaction introduces a competition or 
frustration. The T — k phase diagram for this model shows a rich behavior. It is 
easy to prove analytically that [19j [20l [21] [22] at zero temperature, the system is in 

a ferromagnetic state for k < 0.5, and in antiphase (+ H h H ••■ along x 

direction and all like spins along y direction) for k > 0.5 with a "multiphase" state at 
k = 0.5. (The multiphase state comprises of all possible configurations that have no 
domain of length 1 along the x direction.) But the finite temperature phase diagram 
is yet to be solved analytically. There had been many predictions from various kinds 
of numerical calculations. But here we present the qualitative phase diagram (Fig. [1]) 
that has been obtained by extensive Monte Carlo studies [23l [24] . As the previously 
conjectured floating phase was not found in 23, 24|j it has not been shown in Fig. [I] 

The phase diagram consists of ferromagnetic phase (for k < 0.5) and antiphase 
(for k > 0.5) at low temperature along with a paramagnetic phase at high temperature 
(for all k values). So for both k < 0.5 and k > 0.5 there is a transition temperature 
depending upon the value of k. For k < 0.5, the study of percolation transition is 
done in a similar way as that for the unfrustrated 2d Ising model. But for k > 0.5, 
the cluster cannot be defined as we have for k < 0.5. Because for k > 0.5 below the 
static transition temperature we mainly have striped states with like spins oriented 
along V-axis and a periodic modulation along the direction of frustration. Here we 
propose an alternative definition of cluster to study the percolation transition. By 
a cluster we mean a set of nearest neighbour spins which is in the same state as it 
was in its initial condition. Starting from an exactly antiphase state (initial cluster 
size is the system size) we acheive a steady state at a particular temperature (due 
to thermal agitation some of the spins flip and thus patches of clusters form). Now 
we apply a field pulse which varies spatially along the system such that the field acts 
on the spin along the direction opposite to that with which it started the dynamics 
i.e. sgn(h p ) = sgn(s x . y (0)), where h p is the applied field strength at site x,y and 
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s x,y(Q) is the initial spin value at lattice site x, y. This means that if we start with 

a "+ H h H • ■ ■" spin configuration along cc direction, the applied field will 

be along downward direction for the first two columns and upward along the next 
two columns and this pattern continues over the entire system. As soon as we switch 
on the field pulse (applied along the reverse direction of the initial spin orientation), 
the spins will tend to orient along the field direction and the spins which are at a 
state identical to that of its initial state may flip (depending on the field strength and 
duration) and thus the cluster sizes eventually decreases giving rise to a percolation 
transition. 

In this paper we study the percolation transition for a wide range of frustration 
parameter (0.0 < K < 1.0) and find how the critical field strength varies with the 
frustration parameter for a particular pulse duration. We have used Monte Carlo 
simulation to obtain the numerical results and fully periodic boundary condition 
throught our entire study. We have calculated the critical exponents numerically by 
finite size scaling analysis to see whether the exponents are affected by the frustration 
parameter. In Section 2 we discuss the model studied. In Section 3, we present the 
results obtained for the finite size scaling analysis of the percolation order parameter 
and in Section 4 the fractal dimension and the scaling analysis of the Binder cumulant. 
In the final section we have discussed some consequences related to our study. The 
main observation of this work is that the critical exponents are independent of the 
frustration parameter and also same as that of the normal 2d Ising model. 



2. The Model 



The model studied here is the 2d ANNNI model under a time-dependent external 
magnetic field pulse, described by the Hamiltonian 

H = — J 8 x,y[8te+l,y + S x ,y+1 — KSx+2,y] — h(t) S x>y (2) 
x,y x,y 

where the sites (x,y) runs over a square lattice, the spins s x , y are ±1, J is the nearest- 
neighbor interaction strength (we keep J = 1.0 here on) and n is the competition or 
frustration parameter. The static critical temperature (T®) depends on the value of 
k (Fig. [TJ. For k = 0, T c ° = 2/ln(l + y/2) 2.269.... For k < 0.5, the field pulse 
applied is spatially uniform (directed opposite to that of the initial spin configuration 
we start with) having a time dependence as : 

h{t) = { -\ .*o <*_<*> + At 
w [ otherwise v ' 

For k > 0.5, the field pulse applied is such that the field direction is opposite to that 
of the starting configuration (antiphase order) of the spins, but the time dependence 
is the same as that of Eq. [31 In our simulation we ensure that the time to at which the 
pulse is applied, the system reaches its equilibrium configuration at that temperature 
T(<T C ). 



The spins are selected randomly and are flipped following the normal Metropolis 
algorithm. The energy difference (AE 1 ) to flip the spin is calculated and flipped with 
probability min{l, exp(— AE/T)}. At t = to, the system acheives a steady state at the 
corresponding temperature (in the absence of any field) . As soon as the field pulse is 
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Figure 2. Variations of the largest cluster size (P m ax) with the field pulse 
amplitude (hp) in case of 2d ANNNI model where K = 0.3, T = 0.2 and At = 4. 




Figure 3. The plot of hp against the frustration parameter n. (We have ommitted 
the study close to k = 0.5, because in this region the transition temperatures are 
too low.) 



switched on, a large number of spins which were intact in their initial state, flip during 
the pulse duration. Thus the percolation order parameter P max — Sl/L 2 (where Sl 
is the size of the largest cluster and L is the linear size of the system) decreases during 
this duration and for a particular combination of T, k and At, for a particular value 
of hp, the system undergoes a percolation transition (Fig. [2]). 

It is found that for a particular temperature for k < 0.5, as the value of k is 
increased, the transition field hp decreases, but for k > 0.5, hp increases with k 
(Fig. [3]). For k < 0.5, as n is increased, the increasing frustration itself destabilises the 
order and thus the pulse amplitude needed for transition decreases. But for k > 0.5, 
the situation is different. The increasing frustration parameter stabilises the antiphase 
more and more, thus requiring higher pulse amplitude for the percolation transition 
(Fig.E). 




Figure 4. Inset : The scaled largest cluster sizG i^Pmax LP/ V ) has been plotted 
with h p for six different system sizes (L = 100, 200, 300, 400, 800 and 1000) for a 
typical value of re < 0.5. The crossing point determines the critical field amplitude 
(hp). Main : The value of hp and has been used to make the plots for the 
scaled variables PmaxL^ / " and (h p — h v )L?-l 1 ' collapse by tuning the value of X/v 
for all the different L values. 




Figure 5. Inset : The scaled largest cluster size (PmaxL^ I") has been plotted 
with h p for six different system sizes (L = 100, 200, 300, 400, 800 and 1000) for a 
typical value of re > 0.5. The crossing point determines the critical field amplitude 
(hp). Main : The value of hp and fi/v has been used to make the plots for the 
scaled variables P m axL^ '> " and (hp — hp)L}-l v collapse by tuning the value of X/v 
for all the different L values.. 



3. Finite size scaling analysis of the percolation order parameter 

In this section we discuss about the percolation transition behaviour. For a particular 
combination of k, T and At, we study the variation of P ma x with h p . The percolation 
transition is characterised by several exponents which also determine the universality 
class. The order parameter i.e., the relative size of the largest cluster varies as, 

Pmax~(K-h p f. (4) 
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Figure 6. The plot of hp against At for three different combinations of (re,T) 
for the regime n < 0.5. 
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Figure 7. The plot of hp against At for three different combinations of (re, T) 
for the regime n > 0.5. 



(5) 



The correlation length diverges near the percolation transition point as 

t~( h ;-h p )- v , 

where, h c p is the critical field amplitude. The other exponents can be calculated from 
various scaling relations [25]. But due to finite-size effects it is difficult to determine 
these exponents precisely. We have determined these exponents from the finite size 
scaling relations. P max follows the scaling form 



Pn 



L 1 '" {h c p - h p ) 



(6) 



where J 7 is a suitable scaling function. By tuning fl/v, all the P max L^ ' h '-h p curves can 
be made to cross at a single point and this point determines the critical field amplitude 
(hp). Now once fi/v and h p has been estimated, by tuning 1/u one can collapse 
all the P max L^l u -(h c v — h p )L 1 / u curves on one another, and thus v is determined. 
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We have shown both these plots for k — 0.2 (k < 0.5) and k — 0.7 (k > 0.5) 
in Fig. |3] and Fig. [5] respectively. For both the curves we get the same values for 
P/v = 0.20 ± 0.05 and \ jv = 0.85 ± 0.05 and more importantly this value is the same 
as that we obtained for the normal 2d Ising model found earlier [T8] which indicates 
that the dynamical percolation transition of the 2d ANNNI model also belongs to the 
same universality class as that of the normal 2d Ising model. Only the critical field 
amplitude changes (as shown in Fig. [3]). We have repeated the same exercise for some 
other combination of k, T and At, but the exponents remain unchanged. This suggests 
that the magnitude of frustration does not have any drastic effect on the dynamical 
percolation phenomenon except the quantitative shift of the transition point. 

To have a picture of the effect of various parameters on the critical field amplitude 
we have given two graphs, one for k < 0.5 (Fig. and the other for k > 0.5 (Fig. [7]), 
each showing the variation of hp with At for different sets of k and T . It is quite clear 
from Fig. Othat, for any combination of k and T, h p decreases with At. Also for any 
comination of k and At, hp decreases with T, because thermal excitation always helps 
in percolation transition and thus reducing the necessary external field amplitude. But 
the frustration parameter plays an opposite role on the two halves of the k regime. 
For a fixed pair of T and At, when we are below k < 0.5, with increase of k we need 
less field for reaching the transition point because increasing the frustration parameter 
itself reduces the parallel alignment between the adjacent spins and thus facilitating 
the transition. But above n — 0.5, the increase in frustration parameter supports the 
antiphase structure (with which we begin the dynamics for n > 0.5), and thus we 
have to apply more field amplitude to flip the spins to opposite direction. Thus the 
frustration parameter plays opposite role in cluster formation in the two regimes. 



4. Fractal dimension and Binder cumulant 



At the transition point, the percolation clusters become self similar. The largest cluster 
depends on the linear system size and it varies as Sl ~ L D where Sl is the largest 
cluster and D is the fractal dimension. The fractal dimension measured (Fig. |S]) is 
found to be independent of the value of K and the value is D = 1.82 ± 0.01 which is 
again the same as that we found in case of the normal 2d Ising model. This fractal 
dimension is known to be related to the spatial dimension (d) as D = d — ft] v. This 
relation is also satisfied for our estimated values of f3/v = 0.20 ± 0.05 (mentioned in 
previous section) and D for all values of k. 

For further verification of the critical points and the exponents, we have also 
studied the Binder cumulant of the order parameter, defined as [26] . 

(P 4 ) 

t j -i \ max/ 

U ~ 1 Q/p2 \2 ' 
°\ max/ 

where P m ax is the percolation order parameter and the angular brackets denote the 
ensemble average. The crossing point of the curves U — h p for different system sizes 
gives the critical field amplitude (hp). The Binder cumulant has also been used to 
determine the correlation length exponent as it follows the scaling form 

U=U((h C p-h p )L^), (8) 
where U is a suitable scaling function. Having obtained hp from the crossing point, if 
we plot U with the scaled field amplitude (h'i — hp)L 1 l v , for a suitable value of v we will 
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Figure 8. Variation of the largest cluster size with the linear system size for two 
different values of re = 0.2 and 0.7 at the critical field for T = 0.2 and At = 4. 
The log-log plot gives the fractal dimension to be D = 1.82 ± 0.01 for both the 
values of re. 




Figure 9. Main : The Binder cumulant has been plotted with h p for six different 
system sizes (L = 200, 300, 400, 800 and 1000) for a typical value of re < 0.5. The 
crossing point determines the critical field amplitude (hp). Inset : The value of 
has been used to make the plots for U and (h^ — /ip)!/ 1 /" collapse by tuning 
the value of 1/v for all the different L values. 



obtain a data collapse and that value of v determines the correlation length exponent. 
We have studied these results both for k < 0.5 (Fig. and k > 0.5 (Fig. ITU1). 

For all values of k (only two sets of K, T has been presented in Fig. O and Fig. [TU)) , 
the value of 1/v remains same (0.85 ± 0.05) and also equal to that we have already 
obtained from that of the finite size scaling of percolation order parameter. The value 
of the Binder cumulant at the crossing point for various system sizes is also identical 
(0.62 ± 0.01). The only obvious quantitative change is seen for the transition point. 
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Figure 10. Main : The Binder cumulant has been plotted with hp for six different 
system sizes (L = 200, 300, 400, 800 and 1000) for a typical value of k > 0.5. The 
crossing point determines the critical field amplitude (hp). Inset : The obtained 
value of hp has been used to make the plots for U and (hp — hp)L^'" collapse by 
tuning the value of \jv for all the different L values. 



5. Discussions 

Spatially modulated phases of metallic alloys can be investigated in terms of 
ANNNI model and considering the application of these alloys in magnetisation 
reversal phenomenon, we have studied the magnetic field-pulse induced percolation 
phenomenon in 2d ANNNI model. We have focused on the role of the frustration 
parameter in affecting the percolation transition of the 2d ANNNI model. By applying 
magnetic field-pulse of finite duration (At) for a certain frustration parameter (k) and 
temperature (T) (below the corresponding static transition temperature), we have 
determined the critical field strength (hp) for which the transition occurs. Although 
we have found that the frustration parameter highly enriches the static phase diagram 
of 2d ANNNI model (various novel phases and transitions appear due to k), quite 
surprisingly it does not have any large impact on the dynamical percolation transition 
apart from shifting the transition points. We have restricted our study in the regime 
0.0 < k < 1.0. From finite size scaling analysis we have showed that both the 
percolation order parameter and the Binder cumulant leads to the same values of 
critical exponents as we have found for the normal 2d Ising model (thus it belongs 
to the same universality class as that of the dynamical percolation transition of 
2d Ising model). The values of exponents determined are fi/v = 0.20 ± 0.05 and 
1/v = 0.85 ± 0.05. The fractal dimension D = 1.82 ± 0.01 also confirms the relations 
between the exponents. 

Although we have not presented the details, in this context we have also 
investigated the transition behaviour of the magnetisation (for k < 0.5) and sub-lattice 
magnetisation (for n > 0.5). The sub-lattice magnetisation is defined as, x-direction, 

i — — l 

the sub-lattice magnetization is defined by m s = ~ X^=o ( m 4x+i + W4a:+2 — "mAx+z — 
m 42+4) with, m x = -£ Yl y =i s x,y where L is the length of the system. It has been 
found that the transition points were same as that of the percolation transition (for 
a certain set of value of k,T and At). The study of finite-size scaling of the order 
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parameters also gives the same critical exponents {fijv and l/v). Only the value 
of the crossing point of the Binder cumulant for magnetisation shows non-universal 
behaviour for different temperatures. This value (U* = 0.62 ± 0.01) is robust in case 
of percolation transition. 

This study may be useful in the application of magnetisation reversal of spatially 
modulated magnetic materials where the only change necessary is the value of the 
magnitude of the critical field, keeping the nature of the transition intact. Although 
we have limited our regime of investigation below k — 1.0, one can also explore this 
transition beyond k = 1.0 or other modulated structures. 
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